A microscopic model for spiral ordering along (110) on the MnSi lattice 
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We study an extended Heisenberg model on the MnSi lattice. In the cubic B20 crystal structure 
of MnSi, Mn atoms form lattices of of corner-shared equilateral triangles. We find an ubiquitous 
spiral ordering along (110) for Ji < 0, and J2 = Js > 0, where J\,Ji, and J3 are 1st, 2nd and 
3rd nearest neighbor Heisenberg interactions, respectively. While the ordering direction of (110) is 
reasonably robust to the presence of the Dzyaloshinskii-Moriya interaction, it can be shifted to the 
(111) direction with the introduction of a magnetic anisotropy term for small Jil\J\\. We discuss 
the possible relevance of these results to the partially ordered state recently reported in MnSi. 

PACS numbers: 75.10.Hk,75.25.+z,71.10.Hf 



Introduction: The helical magnetic ordering in MnSi 
at ambient pressures has been understood as a result of 
Dzyaloshinskii-Moriya (DM) interactions which arise due 
to the lack of inversion symmetry of the lattice 0. Pre- 
vious studies of symmetry respecting Landau-Ginzburg 
(LG) free energies have successfully described phenom- 
ena observed in 

MnSi[lH00 

including the wavevector 
reorientations induced by a magnetic field 0]. In the ab- 
sence of an applied field, LG theories showed that ferro- 
magnetic (FM) and DM terms lead to a helical ordering 
with a fixed periodicity, but no orientational preference. 
The spiral ordering along (111) seen at ambient pressures 
can been understood as arising from crystal anisotropy 
terms consistent with the lattice symmetry which 
give energetic minima either at (111) or (100). However, 
ordering along (110) cannot be obtained from an analysis 
of LG free energies solely consisting of the FM, DM, and 
crystal anisotropy terms usually considered in MnSi. 

It was thus of great interest when neutron scattering^ 
on MnSi under high pressures revealed unusual magnetic 
properties. At pressures above p c = 14.6kbar, the neu- 
tron scattering intensity no longer sharply peaked along 
(111). Instead, broad peaks at (110) with more diffuse 
scattering over a sphere of fixed magnitude wavevector 
have been reported in a small window of the pressure 
and temperature phase diagram. The regime where this 
phenomenon has been observed has been dubbed the par- 
tially ordered state of MnSi. The magnitude of the order- 
ing wavevector has not dramatically changed across p c , 
but its orientation has clearly shifted from (111) to pre- 
dominantly (110). It also remains to be understood why 
the resistivity has been found0,E| to be proportional to 
T2 at and well above p c , which was not anticipated by 
itinerant models Q. 

Recently, to describe the partially ordered state, 
the addition of a pseudoscalar order parameter was 
argued^ to lead to a "blue quantum fog" with proper- 
ties of a chiral liquid. Attempts to explain the observed 
shift of weight from the (111) direction of the ambient 
pressure ordered helimagnet to the (110) direction at high 
pressures has led to theories which: i) favor a weakly or- 
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FIG. 1: The 3D corner-shared equilateral triangle 
ftrillium lTlln lattice, (a) Nearest neighbor (nn) bonds (Ji) 
form equilateral triangles with normal vectors along relative 
(1,1,1) directions, C3 symmetry axes of the structure. Each 
point (labeled a, j3, 7, 8) belongs to three triangles. The sub- 
script i denotes the elements of the cubic unit cell as given 
before Eq0 (b) Second and (c) third nn bonds again form 
trillium lattices with normal vectors along the relative (1,1,1) 
directions, (d) Looking down the (1,1,1) axis, the 1st, 2nd 
and 3rd nn bonds and their relative distances and orienta- 
tions are shown, "x" marks the location of the nn 8 sites one 
layer down and "o" one layer up. Although the ground state 
of AF couplings on a trillium lattice features 120° rotated 
spins, it is not possible to simultaneously satisfy any two of 
{Ji, J2, J3}. 



dered state resulting from introducing (non-local) mode- 
mode coupling interactions between helimagnetic struc- 
tures to generate a new type of spin crystal [T3| : or ii) 
favor a lattice of cylindrical skyrmionic tubes with axes 
pinned to the (111) direction|13|. 

However, the similarities between the MnSi lat- 
tice structure (B20) and the corner-shared tetrahcdral 
(L1V2O4, (Yo.97Sco.o3)Mn 2 ) and distorted windmill (/?- 
Mn) lattices have been overlooked so far. Heisenberg 
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FIG. 2: Dispersion curves for J2/IJ1I = J3/I Ji = 1/3 and |Ji| + J2 + J3 = 1 along the (100), (110) and (111) directions. For 
all curves there is a mirror symmetry (not shown) about the left and right axes and we set ao=l. The ground state has energy 
-2(| Ji I + J2 + J3) with Q = (71", 7T, 0) in agreement with the rotor model. 



models have been used to successfully explain anomalous 
behavior in these itinerant svstems |l4ill5lllq . Moreover, 
the same anomalous resistivity oc Ti near and above the 
pressure tuned quantum critical point has been reported 
111 FeGe[I3, which has the same B20 lattice structure. 
This raises the possibility that the lattice structure may 
play a role in understanding the helical magnetic ordcr- 
ings in MnSi. In this paper, we suggest a microscopic 
model consisting of Heiscnbcrg interactions on the MnSi 
lattice. We find that a spiral ordering along (110) is nat- 
urally realized as a result of a competition between FM 
and antiferromagnetic (AFM) interactions. 

Model: The unit cell of MnSi is a simple cubic struc- 
ture, and contains four atoms at sites: {Si, on, Pi} = 
{(u, u, u); (u + |, I — it, 1 — it); (1 — it, u + |, I — u); (k — 
u, 1 — it, )}ao with UMn = 0.138 and ao = 4.56A|l7||. 
Si atoms are located in the same pattern with ust = 
0.854. Mn atoms form corner-sharing equilateral trian- 
gular (hereafter trillium[llT|) lattices as shown in Fig. ^ 



Each Mn atom is shared by 3 triangles and the normal 
vectors to the triangles lie along the 4 body diagonals of 
the cubic unit cell: (1,1,1), (1,-1,-1), (-1,-1,1) and (-1,1,- 
1). We study a classical Heiscnbcrg model with near- 
est neighbor (nn) FM interactions (Ji < 0) and next- 
( J2 > 0) and third- ( J 3 > 0) nn AFM couplings on the 3 
dimensional (3D) trillium lattices formed by the Mn sites 
of MnSi, 

H = Ji ~~Sj ■ ~Sj + Ji ~s t ■ ~~Sj + J3 2J ' 

(«> <(«•» mm 

C=i C=i Q 

Here represents the £th order Heisenberg coupling 
with £=1-3 and S q — C^qt's'q, ~^q) where ~Si — 
J2q ~s V Qe % Q Ti . Mq is the matrix given by, 



/ cos(-qq x a )e l ^ q y a ° +vq * a ^ 

cos(77 fe a )e- l (^ ao+ ^^ a °) 

cos{m y ao)e- l{u ^ ao+ ^ ao) cos(i 1 q z a )e- 1 ^ a °~' /q y a ^ 

\ cos(r]q z ao)e- i( -^ ao+l/q y a ^ cos^ao^-^ 00 -^^ 



cos(r 1 q y ao)e i ^ q * aa+ w* ao *> cos(m z ao)e i( -^ a ° +yq y a ^ \ 

cos(f]q z a )e l( -^ a «- l ' q y a «'> cos(r]q y ao)e t( - l/ ^ ac >-f iq ' ao '> 

cos(7 1 q x a )e l ^ q y a °-' yq ' ao) 
cos(rjq x a )e- l ^ q y a °- vq ' a ^ / 



where 77= ( 2 "+i) ; ^=(2u+m— i), v={2u+p) and the three 
integers {m, n,p} take the values {0, 0, 0}, {0, 0, —1} and 
{1, 0, 0} for ( = 1, 2 and 3 respectively. 

A mean field (MF) treatment which corresponds to 
the diagonalization of the matrix of Mq^ has been 
studied [l9|. The ground state of with energy -3Ji 
shows degenerate states satisfying the equation, 

9 

cos 2 (rjq x a ) + cos 2 (rjq y a ) + cos 2 (r]q z a ) = - (3) 



Here Q{= {q x , q y , q z )) is a wavevector of the first Bril- 
louin zone and the MF ground state is given by the mini- 
mum of the eigenvalues Xq at E = Aq q where £ represents 
4 different eigenvalues, and Qo is where the minimum oc- 
curs in Q. The set of wavevectors satisfying Eq. Inform a 
sphere-like shape which resembles the sphere of wavevec- 
tors observed in MnSi, while the peculiar anisotropic in- 
tensities cannot be captured. However, a Monte Carlo 
treatment |2dl| on this lattice found an ordered state in 
agreement with the result obtained by a classical rigid ro- 



3 



tor model^3- The ground state of forms a lattice of 
120° rotated spins on each triangle with the wavevector 
along the (100) direction. It is interesting to note that, 
unlike on the trillium lattice, MF and Monte Carlo results 
agree on the corner-shared tetrahedral lattice, common 
to the pyrochlore and spinel structures.. 

It is straightforward to see that the solely 2nd (or 
3rd) nn Heisenberg model generically has the same 
structure|2l| and ground state as the 1st nn model. How- 
ever, as indicated in Fig. ^ it is not possible to simul- 
taneously satisfy AFM couplings for any two of these 
lattices. Below we will show that the microscopic model 
of Eq. ^ finds a ground state with spiral ordering along 
(110) due to the competition between AFM couplings 
and between AFM and FM couplings. We use both MF 
and the rigid rotor approximations to avoid any artificial 
effects of the MF approximation as found in the nearest 
neighbor AFM model. As discussed elsewhere fioj. the 
rigid rotor approximation corresponds to the mapping of 
a spin to a classical vector with 2 independent angles. 
Minimizations of energy with respect to configurations 
of these real space rigid rotors allow one access to the 
physical spin structure of the ground state. 

Phase diagram: We studied the energy dispersion 
curves in Q within the MF approximation for a given 
J 2 /\Ji\, and we set J 2 = J3 for all cases. We found the 
ground state to be FM for J 2 /\Ji\ < 0.214. As J 2 (= J 3 ) 
increases, spiral ordering with a finite wave vector be- 
comes the ground state. The energy dispersion curves 
for three different directions, (100), (110), and (111), are 
shown in Fig. E|for Ji/\Ji\ = 1/3. As shown, we found 
a ground state with energy -2(| Ji| + J 2 + J3) along the 
(110) direction with magnitude Q\ = \ Q | of V2ir/ao rad. 
The rotor model minimizes the energy along the (110) di- 
rection for the same value of J 2 /| Ji| (=1/3) verifying the 
MF result. We present the ground state energies of the 
rotor model as a function of J2/ 1 Ji | for three different di- 
rections in Fig. 13(a). It is clearly shown that magnetic 
ordering along (110) has a lower energy than along the 
(100) and (111) directions for virtually the entire phase 
diagram with exceptions close to Jf"(= 0.214| Ji|) and 
J2/ \Ji I — > 00. Close to Jf", various states represented 
by different directions of the ordering wavevectors are 
almost degenerate. We also present the magnitude of 
the ordering wavevector Qi = |Q| for the corresponding 
states in Fig. [3] (b). The MF and rotor model approxi- 
mations are found to agree for 0.32< J 2 /|Ji| < 0.5. How- 
ever, the MF approximation overestimates the value of 
Qi compared to the rotor model for Jf™ < J 2 < 0.32 1 Ji|, 
while it underestimates the value of Qi for 0.5| J\\ < J 2 . 
It is also interesting to note that the size of Qi for the 
(110) direction is generally larger than those of the other 
two directions. 

Spin structures: As we discussed, an advantage of the 
rotor model is that we can identify the spin structure 




FIG. 3: For Ji < 0, J 2 = J s > 0, and | J r \ + J 2 + Jz = 1 and 
three different directions (100), (110) and (111): (a) Energy 
(in units of (| Ji\+J 2 +Js)) versus J 2 /\Ji\. For J 2 /\Ji\ < 0.214 
the ground state is FM with energy 6(Ji + J 2 + J3). For 
J 2 /\Ji\ > 0.214, the ground state orders along (110), except 
for tiny regions of (100) ordering very close to FM (see in- 
set), and (111) order when J 2 /\J±\ —* 00. (b) Magnitude 
of the ordering wavevector, Qi, at the minimal energy solu- 
tion of a given direction versus J 2 /\Ji\. Notice that the min- 
ima have similar Qi magnitudes near the onset of FM. Mean 
field theory tends to over/ (under) estimate Qi for FM/(AFM) 
dominated couplings for the (110) and (111) symmetries. For 
the (110) ordered ground state, a small component of the 
spins (orthogonal to the component rotating with Qi) at the 
{a,7}/({<5, 0\) sites retains a (w/ao, ir/ao, 0) symmetry. 



of a given state. Spin structures periodic along (110) 
tend to break the four-site unit cell into two sublattices. 
Indeed, all ground state spin structures with this sym- 
metry feature parallel spins on the nn /3, and Si+y sites 
(here y denotes a coordinate increase by one unit in the 
y-direction). To illustrate this, let us consider the special 
point, J 1 = — J 2 — J3, J 2 = J3 > 0, where a further choice 
of parallel spins at the ji and sites minimizes the en- 
ergy at — 2(1 Ji| + J 2 + J3) with wavevector ( — ,— ,0). 
In this configuration, two spins on each FM (Ji) bonded 
triangle are parallel, while two spins on each AFM (J 2 
and J3) bonded triangle are anti-parallel as shown Fig. 
0| The energy is independent of the relative orientation 
of the third spin on each triangle, as FM and AFM inter- 
actions cancel, leaving the structure independent of an 
overall relative angle between pairs of spins. 
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FIG. 4: Looking down the (111) axis. At the special point 
Ji = — 2J2 = — 2J3, the ground state spin structure has (110) 
symmetry and spin arrangements form two independent sets 
of pairs {a, 7} and {S, /3} (ie. the relative angle, 9, between 
these two sets is completely arbitrary.) Two spins on each 
FM (Ji) bonded triangle are parallel, while two spins on each 
AFM (J2 and J3) bonded triangle are anti-parallel. An anal- 
ogous state is not available to (100) and (111) symmetry spin 
structures. 



For a coplanar spiral spin structure (even one with 
structure inside the unit cell as here), all spin compo- 
nents should lie in a plane which, for convenience, here 

— y 

we will define to be perpendicular to the wavevector Q 1 . 
This is the case for 0.32 < J 2 /\Ji\ < 0.5 where the 
MF and rotor models agree. For a 3D spin structure, 
a non-coplanar spin structure may be energetically fa- 
vorable, where one or two spins of the four site unit cell 
slightly tilt out of this plane. Further, the spin compo- 
nent out of the plane can alternate in sign, exhibiting 
a distinct wavevector, Qi- The two sublattice nature 
of the spin structure continues away from the special 
point at J 2 /|Ji|(= J 3 /\Ji\) = 1/2. For J 2 /\.h\ < 0.32, 
spins at the a and 7 sites acquire a component paral- 
lel to Qi with wavevector Q2 = {ft/aa, 7r/ao, 0) ^ Q±. 
For J2/IJ1I > 0.5 a similar noncoplanar spin structure 
arises but at the (5 and (3 sites. In the calculation of the 
neutron scattering structure factor no interference terms 
arise between the spin components with wavevectors Q 1 
and Q 2 , as the corresponding spin structures are orthog- 
onal. The only effects of the development of this second 
ordering wavevector on the neutron scattering structure 
factor would be to slightly decrease the weight at Q 1 , and 
to create a tiny peak at Q2 = (7r/ao, 7r/a , 0). We have 
extracted the components of the spin structure along and 
perpendicular to Qi from the ground state (110) spin 
structures as a function of J2/I J\\ although we omit this 
here for brevity. 

Discussion: It is interesting to ask whether this (110) 
ordered state might have any relevance to the weakly or- 
dered (110) state found by Pflciderer et alQ. We know 
from measurements!^ of magnetic susceptibility, (x), 



that the low temperature Curie law maintains its slope 
as pressure increases past the transition to the anomalous 
NFL phase. The intercept of x _1 appears to eventually 
reach a weakly negative offset consistent with the pres- 
ence of some AFM couplings. Such couplings might arise 
in an RKKY-based picture of itinerant magnetism be- 
tween local moments on the Mn sites. The pitch of the 
structured helimagnet found here is an order of magni- 
tude larger than that observed in MnSi. This can be 
explained by the inclusion of further neighbor couplings 
which also form trillium lattices. We expect that in this 
way an analogous ordering along (110) can be obtained 
with a smaller wavevector, although the derivation of the 
strength and size of such couplings remains to be inves- 
tigated. 

At ambient pressure, a good approximation to the 
ground state of MnSi is given by a simple spin spiral 
with (111) symmetry. This symmetry pins at low en- 
ergies from a nearly degenerate set of fixed period sim- 
ple spirals sampled at higher temperatures. The lack 
of inversion symmetry of the MnSi lattice permits the 
introduction of asymmetric spin terms arising from the 
spin-orbit coupling on the lattice. The addition of DM- 
vectors parallel to the nearest neighbor Mn-Mn bonds on 
this lattice nicely reproduces a nearly degenerate surface 
in momentum space as will be presented elsewhere |23j. 
With the noted exception of the FM phase (which be- 
comes a (100) ordered phase), we find that the phase 
diagram of Fig. |3|is very weakly effected by competition 
from these first-neighbor DM-vectors. 

Briefly we comment on how the (111) ordered ambi- 
ent pressure phase is thought to pin from the degener- 
ate helimagnetic phase. The local environment of the 
Mn site is surrounded by 7 Si atoms, with the shortest 
Mn-Si distance along the local (1,1,1) axis. It therefore 
seems plausible that crystal/ligand field effects might 
produce a local (1,1,1) anisotropy of the spin structure. 
One term which might arise from a full consideration of 



spin-orbit interactions is of the form, — £ J2i v ( s 



i : Ising 



Here the sum runs over the number of unit cells and the 



4 sites within the unit cell, s 



i. I sing 



( s i ■ ~e u )~e u and 



e v is the unit vector associated with the local (1,1,1) 
anisotropy at site v. Assuming this term is small, we 
can evaluate the energy change that such a term would 
introduce to the simple helimagnet. For ( > this term 
is found to lower the energy of the (100), (110) and (111) 
directions by -6£, -6.5£ and -8£ respectively per unit cell, 
effectively favoring (1,1,1). A pressure increase likely 
increases the overlap of the various Si atoms with the 
Mn sites (as apparently seen in FcGc0]), reducing the 
strength of pinning. 

Summary: In this paper we have provided a simple 
classical spin model on the MnSi lattice which exhibits 
an extended phase with a (110) ordered ground state. We 
have considered an extended Hciscnberg model with both 
FM and AFM couplings joining the Mn sites. Within this 
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microscopic model, we find that the corner-shared equi- 
lateral triangle structure of the lattice plays a crucial 
role in the suppression of both (100) and (111) order- 
ings. Even though AFM couplings are known to favor 
(100) orderj2^|, competitions between AFM J2 and J3 
frustrate this order and are found to yield (110) order 
in the presence of a FM Ji, with or without a nn DM 
interaction. Such physics has not been considered in LG 
studies. A detailed derivation of the RKKY interaction 
on this unusual lattice is desired and could be anticipated 
to lead to a weak (110) ordering arising from a competi- 
tion between AFM and DM interactions. It is tempting 
to speculate that the crossover from (111) to weak (110) 
order seen in MnSi might then result from the application 
of pressure reducing crystal field splittings and increasing 
AFM coupling strengths. 
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